Gravitational wave production: A strong constraint on primordial magnetic fields 
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We compute the gravity waves induced by anisotropic stresses of stochastic primordial magnetic 
fields. The nucleosynthesis bound on gravity waves is then used to derive a limit on the magnetic 
field amplitude as function of the spectral index. The obtained limits are extraordinarily strong: 
If the primordial magnetic field is produced by a causal process, leading to a spectral index n > 2 
on super horizon scales, galactic magnetic fields produced at the electroweak phase transition or 
earlier have to be weaker than B\ < 10~ 27 Gauss! If they are induced during an inflationary phase 
(reheating temperature T ~ 10 15 GeV) with a spectral index n ~ 0, the magnetic field has to be 
weaker than B\ < 10 _39 Gauss! Only very red magnetic field spectra, n ~ —3 are not strongly 
constrained. We also find that a considerable amount of the magnetic field energy is converted into 
gravity waves. 

The gravity wave limit derived in this work rules out most of the proposed processes for primordial 
seeds for the large scale magnetic fields observed in galaxies and clusters. 



PACS Numbers : 98.80.Cq, 98.70.Vc, 98.80.Hw 
I. INTRODUCTION 

Our galaxy, like most other spiral galaxies, is perme- 
ated by a magnetic field of the order of B ~ 10 _6 Gauss. 
Recently, similar magnetic fields have also been observed 
in clusters of galaxies on scales of up to A ~ 0.1 Mpc jjJU. 
There is an ongoing debate whether such fields can be 
produced by charge separation processes during galaxy 
and cluster formation ]3| or whether primordial seed 
fields are needed, which have then been amplified by sim- 
ple adiabatic contraction or by a dynamo mechanism. In 
the first case, seed fields of B ~ 10~ 9 Gauss are needed 
while in the second case B ~ 10 _20 Gauss j| or even 
lCP 30 Gauss in a universe with low mass density Q suf- 
fice. Several mechanisms have been proposed for the ori- 
gin of such seed fields, ranging from inflationary produc- 
tion of magnetic fields J5|-f7| to cosmological phase tran- 
sitions ||. 

Primordial magnetic fields have been constrained in 
the past in various ways mainly by using their effect on 
anisotropics in the cosmic microwave background [p|~[i4| . 
In these works constant magnetic fields and stochastic 
fields with red spectra n ~ — 3 G3] have been consid- 
ered and the limits obtained where of the order of a 
fewx 10 _9 Gauss. A simple order of magnitude estimate 
shows that, from the CMB alone, one cannot expect 
much stronger constraints of magnetic fields: The energy 
density in a magnetic field is 

B 2 

Q B = ^ lO^OvfB/IO^Gauss) 2 , (1) 

8np c 

where i7 7 is the density parameter in photons. We 
naively expect a magnetic field of 10 _8 Gauss to induce 
perturbations in the CMB on the order of I0~ 5 , which 
are just on the level of the observed CMB anisotropies. It 
is thus expected that CMB anisotropies cannot constrain 
primordial magnetic fields to better than a few tenths of 
this amplitude. 



In this work we constrain magnetic fields by the gravity 
waves which they induce classically, via the anisotropic 
stresses in their energy momentum tensor. These grav- 
ity waves lead to much stronger constraints than CMB 
anisotropies, especially for spectral indices n > —3. This 
comes from the fact that the spectrum of the gravity wave 
energy density induced by stochastic magnetic fields is al- 
ways blue (except for n — —3 where it is scale invariant) 
and thus leads to stronger constraints on small scales 
than on the large scales probed by CMB anisotropies. 

The effects of a constant magnetic field on gravity wave 
evolution and production have been studied in ]f6[ | . Here 
we concentrate on the production of gravity waves, but 
consider a stochastic magnetic field. 

The remainder of this paper is organized as follows: 
In the next section we define the initial magnetic field 
spectrum and its evolution in time, and we determine the 
magnetic stress tensor which sources gravity waves. In 
Section 3 we calculate the induced gravity wave spectrum 
and estimate the effect of back-reaction. In Section 4 
we derive limits on the primordial magnetic field using 
the nucleosynthesis limit on gravity waves and discuss 
our conclusions. In order not to loose the flow of the 
arguments, several technical derivations are deferred to 
three appendices. 

We use conformal time which we denote by 77; the scale 
factor is a(rj). Derivatives w.r.t conformal time are de- 
noted by an over-dot, ^ = a. We normalize the scale 
factor today to a(rjo) = 1. The index on a time depen- 
dent variable always indicates today. We assume a spa- 
tially flat universe with vanishing cosmological constant 
throughout. Neglecting a possible cosmological constant 
modifies the evolution of the scale factor only at very late 
times, z < 2 and is therefore irrelevant for the results of 
this paper. We set the speed of light c = I so that times 
and length scales can be given in units of sec, cm or Mpc, 
whatever is convenient. With our conventions, the scale 



f 



factor is given by 



a(ry) = H r)(—^- + v^rad) 



(2) 



where Ho = (3.086 x 10 17 sec) _1 /io is the Hubble param- 
eter, 0.5 < h < 0.8 and f2 ra d = 4.2 x lO -5 ^ 2 is tnc 
radiation density parameter (photons and three types of 
massless neutrinos). 

Note that the scale factor has no units, but conformal 
time and comoving distance do. The normalization of 
a implies that comoving distance becomes physical dis- 
tance today. The conformal time r\ is the comoving size 
of the horizon. The relation between 77 and redshift or 
temperature is simply 



a(r?) 

T{rj) = z(rj)T ~ z(ry)2.4 x 10^ 4 eV 



(3) 



The comoving time of equal matter and radiation, defined 
by a(r) eq )~ 3 = £l ra da(f/ e g)~ 4 or z eq + 1 = Q - ^, is 



rj eq = 2(\/2 - lj^fW^ 1 ~ 1.7 x 10 15 sec. 



(4) 



Greek indices run from to 3, Latin ones from 1 to 3. 
Spatial (3d) vectors are denoted in bold. 



II. PRIMORDIAL STOCHASTIC MAGNETIC 
FIELDS 

In this section we closely follow Ref. During the 
evolution of the universe, the conductivity of the inter- 
galactic medium is effectively infinite. We can decouple 
the time evolution of the magnetic field from its spatial 
structure: B scales like B 2 (rj,x) = B§(x)/a 4 on suffi- 
ciently large scales. (In our coordinate basis Bi on l/a 
and B l oc a~ 3 as can be derived easily from Maxwell's 
equations in curved spacetime with vanishing electric 
field, see e.g. |l7]]). On smaller scales, the interaction 
of the magnetic field with the cosmic plasma becomes 
important, leading mainly to two effects: on interme- 
diate scales, the field oscillates like cos(u Akf]), where 
i'a = (B 2 /(4tt(p + p))) 1 / 2 is the AlfVen velocity, and 
on very small scales, the field is exponentially damped 
due to shear viscosity [|l8-20|. We will take into account 
the time dependent damping scale as a time dependent 
cutoff fed (77) in the spectrum of B. As we shall see, our 
constraints come from small scales where the spectrum 
is exponentially damped and oscillations can be ignored. 
We therefore disregard them in what follows. The ex- 
pressions for kd{rj) are derived in Appendix |A|. The only 
result of this appendix relevant here is that the damping 
scale 1 /kd{rf) grows like a positive power a > of r\ and is 
always smaller than the horizon scale, kd{rf) oc l/r/ a and 
kd{i]) > 1/rj. The reader not interested in the details of 



damping and confident with this relatively obvious result, 
can skip Appendix [A|. 

We model Bo(x) as a statistically homogeneous and 
isotropic random field. The transversal nature of B then 
leads to 

(B l (k)B*i(q)) = <5 3 (k- q)(<f J - k l ti)B 2 {k) . (5) 
We use the Fourier transform conventions 

B j {k)= J d 3 xexp(ix-k)Bl(x) , 

Bg(x) = 7^3 J d 3 fccxp(- lX .k)^(k) , 

and k = k/fe, k = \/J2i(k 1 ) 2 ', k is the wave vector to- 
day which is also the co-moving wave vector. Its unit is 
inverse length which we will express in sec -1 . 

We want to derive a limit on the amplitude of magnetic 
fields on the scale A ~ O.lMpc generated by a primordial 
process which took place before ij — O.lMpc ~ 10 13 sec 
corresponding to T ~ IkeV. Hence we are mainly in- 
terested in magnetic fields generated on super horizon 
scales. As we shall see, our limits only apply for fields 
generated before nucleosynthesis, T > T nuc ~ O.lMeV. 
The main examples we have in mind are inflationary gen- 
eration of magnetic fields (||| j magnetic fields generated 
in string cosmology and magnetic fields generated dur- 
ing the electroweak phase transition || . 

In the first two examples, a simple power law magnetic 
field spectrum with upper cutoff k c ~ 77"^ is generated. 
The conformal time rji n marks the end of inflation or the 
string scale respectively. 

Electroweak magnetic field production is causal, lead- 
ing mainly to fields on scales smaller than the size of 
the horizon at the phase transition, i] ew ~ 4 x 10 4 sec 
~ 10 15 cm ~ 3 x 10 _4 pc. These sub-horizon fields, which 
cannot propagate into larger scales during the linear evo- 
lution discussed in this paper, and which are essentially 
damped by viscosity, will be neglected in this paper. Mo- 
tivated from inflation, we simply impose an initial cut- 
off scale k c (rjin) = l/^m- Allowing for more small scale 
power, as it is certainly present initially in causal mecha- 
nisms, only strengthens our result which actually comes 
from the smallest scales not affected by damping. 

If B is generated by a causal mechanism, it is uncor- 
rected on super horizon scales, 



(B i (x, rj)B j (x',r)))=0 for |x - x' | > 2?y 



(6) 



Here it is important, that the universe is in a stage 
of standard Friedman expansion, so that the comoving 
causal horizon size is about r\. During an inflationary 
phase, the causal horizon diverges and our subsequent 
argument does not apply. In this somewhat misleading 
sense, one calls inflationary perturbations 'a-causal'. 

According to Eq. (@), (B i (x,ri)B j (x' ,r?)> is a function 
with compact support and hence its Fourier transform is 
analytic. The function 



2 



(£ l (k)B* J (k)) = (5 l J - k l k j )B 2 {k) 



(7) 



is analytic in k. If we assume also that B 2 (k) can be 
approximated by a simple power law, we must conclude 
that B 2 (k) cx k n , where n > 2 is a even integer. (A 
white noise spectrum, n — does not work because of 
the transversality condition which has led to the non- 
analytic pre- factor 5^ — k 1 ^.) By causality, there can 
be no deviations from this law on scales larger than the 
horizon size at formation, r\i n . As explained above, we 
neglect fields on smaller scales by a simple cutoff. 

We assume that Bo is a Gaussian random field. Al- 
though this is not the most general case, it greatly sim- 
plifies calculations and gives us a good idea of what to 
expect in more general situations. 

Using Wick's theorem for Gaussian fields we can cal- 
culate the correlator of the tensor contribution to the 
anisotropic stresses induced by the magnetic field, which 
we denote by Ely . One finds (see Appendix pi) 



(ir j (k,7 7 )ir im (k',77)) = \Ii(k,7 1 )\ 2 /a 12 M l3lm 5{\i-\ ; !) 

4 

, . . M,- ..ill ■■■'.> i ],•' , , \ \ — 



n^k^ir^k',^)) = -g/ 2 (MWk-k'), (8) 



where 

M ijlm (k) ee 6 il 8 jm + 8" n 8 ]l - S lj S lm + k- 2 (5 l] k l k m + 
S lm k l k j - S'ttfk" 1 - S" n k l k j - S jl k l k m 
~5 ]m k l k l ) + k~ A k i yk l k m , (9) 

and 

/(fc)2 = 16(2^ / d3qB2{q)B2 ^ k - q|)(1 + 2 ^ + ^ ' 

(10) 

with 7 = k • q and [3 = k • k — q. For this result we 
made use of statistical isotropy, which implies that the 
two spin degrees of freedom of IEy have the same average 
amplitude. More explicitly: in a coordinate system where 
k is parallel to the z-axis, Ely has the form 



together with Eq. (g) , statistical isotropy then gives 




|n + | 2 )H|n x | 2 ) 



r/ 2 



(11) 



To continue, we have to specify B 2 (k). For simplicity 
we assume a simple power law with cutoff k c which can 
depend on time. As all scales smaller than \/kd{r/) are 
damped, clearly we have to require k c (rf) < kd(rf). Moti- 
vated by inflationary magnetic field production we choose 
k c (Vin) ~ l/?7m, the primordial magnetic field is coherent 
up to the horizon size at formation. For magnetic fields 



produced during the electroweak phase transition, the 
'coherence scale' is substantially smaller |21 , k c (rji n ) ^> 
l/f]i n which would strengthen our limit as we shall see. 
Since it is unphysical to assume k c (rj in ) < l/r}i n , our 
assumption is conservative. We set 

k c (r]) = mva(l/r)i n ,kd('r])) ■ 

It is important to keep in mind, that this cutoff scale is 
always smaller than the horizon scale. 
We now can parameterize B 2 by 

f (2tt) 5 (A/V2)"+ 3 p2,n t nr k<-k 

B 2 (k) = \ 2 r[nj2] X lorfc<fc c (12) 
[ otherwise. 

The normalization is such that 

B l = yJ d 3 r(S (x)i? (x + r))exp(-^) , (13) 

where V = J d 3 r exp(— r 2 /2A 2 ) = A 3 (27r) 3 / 2 is the nor- 
malization volume. (We have assumed that the cut- 
off scale is smaller than A.) We will finally fix A = 
0.1/i _1 Mpc, the largest scale on which coherent magnetic 
fields have been observed; but the scaling of our results 
with A will remain obvious. 

The energy density in the magnetic field at some ar- 
bitrary scale £ is cx B 2 cx B 2 (k)k 3 \ k=1/e cx £~( n + 3 ). In 
order not to over-produce long range coherent fields, we 
must require n > —3. For n = — 3 we obtain a scale 
invariant magnetic field energy spectrum. 

Using Eqs. (|l2| ) and (|o|) we can calculate /. The inte- 
gral cannot be computed analytically, but the following 
result is a good approximation for all wave numbers k E3] 



f 2 (k, V )~Ax 



with 



fc c (r?) 2 ™+ 3 
fc 2n+3 



for n > -3/2 
for n < -3/2 



(14) 



A = (2tt) 3 (\/V2) 2n+e Bi 



1G 



For n > —3/2, the gravity wave source II is white noise, 
independent of k. Only the amplitude, which is propor- 
tional to (Afc c ) 2n , depends on the spectral index. This 
is due to the fact that the integral (|l0|) is dominated by 
the contribution from the smallest scale k^ 1 . The in- 
duced gravity wave spectrum will therefore be a white 
noise spectrum for all n > —3/2. 



III. GRAVITY WAVES FROM MAGNETIC 
FIELDS 



We now proceed to calculate the gravity waves induced 
by the magnetic field stress tensor. The metric element 
of the perturbed Friedman universe is given by 

ds 2 = a 2 (r))[dr] 2 - (<% + 2h i3 )dx l dx 3 } , 



3 



where h\ = and h\k % =0 for tensor perturbations. The 
magnetic field sources the evolution of hij through 



ha + 2-hi 



k 2 hij = 8ivGIlij 



(15) 



Hij is a random variable, but its time evolution is deter- 
ministic, it evolves in time simply by redshifting and by 
the evolution of the cutoff. Each component is given by 

n.(M) = 4/(M)n.(*0 , 

where II, (k) is a time independent random variable with 
power spectrum (|II.(fc)| 2 ) = 1. Therefore, also each 
component of the induced gravity wave is given by 

h.{k,ri) = h(k,r])U.(k) , 

where h(k, rf) is a solution of 

8ttG 



h + 2-h + k 2 h= ^-f{k lV ) 
a °> (.V) 



(16) 



The gravity wave power spectrum is then given by 

(fe«(k, f? )^.(k',» 7 )) =Ah 2 (k,r,)S(k-k') . (17) 
In real space ,the energy density in gravity waves is 



PG 



IGirGa 2 ' 



The factor 1/a 2 comes from the fact that h denotes the 
derivative w.r.t. conformal time. Fourier transforming 
this relation, we obtain with Eq. ( |l7| ) 



(>G 



with 



K dk dp G (k) 
k dlog(fc) 

dp G (k) 



(18) 



k 3 h 2 



such that 

dtt G (k) 



dlog(fc) a 2 (27r) 6 G 
dp G (k) 



k 3 h 2 



dlog(fc) p c d\og{k) a 2 p c {2Tr) e G 



(19) 



where p c = 3H 2 / (8nG) denotes the critical density to- 
day. In Appendix ^| we solve Eq. ( |l6| ) for n < —3/2, when 
/ is time independent, and we show that for wave num- 
bers which enter the horizon in the radiation dominated 
era, the density parameter in gravity waves produced by 
the magnetic field can be expressed as 



dfl G (k) 12fc 3 f(k) 2 log 2 (x in ) 
dlog(fc) ~ p 2 O rad (2^)5 
for n < -3/2 . 



(20) 



Fourier transforming the expression for the magnetic field 
energy p B = (B 2 (x)) / (8tt) , we obtain the magnetic field 
density parameter at time rj, 



<m B {k) 
fi B (fc c fa)) 



(k\) 



n+3 



8n Pc 2(»+3)/2r(^±3) 

k " {7l) dk dn B {k) 



k dlog(fc) 
(k c X) n + 3 



8tt Pc 2(«+5)/2r(H±5) 



(21) 



(22) 



Note that £l B may well be considerable on small scales, 
since this is the magnetic field energy at very early times 
which can be damped and transformed, e.g. into radi- 
ation later. But of course, for our perturbative calcu- 
lation to apply, we must require ^,g(fc) < fl ra d during 
the radiation dominated era. Using Eqs. j2l|j22] ) and the 
result (14|) for /, we obtain from Eq. (|2C 



dQ G (k) 
d^gjk) 



241og 2 (x m ) , 



n G = 



^rad 

for - 3 < n < -3/2 

lhm dk dn G (k) 

o 

2 B y 



(23) 



Q B (j]in) 



k dlog(fc) 
12(ra + 3) 



^rad 

for - 3 < n < -3/2 



(24) 



In the integrated formula for fl G we have neglected the 
logarithmic dependence \og 2 (xi n ). 

If n > —3/2 the result changes since / now depends 
on time via the cutoff k c {rf) = min(l/?7i n , kd(i]))- Clearly, 
kd(ijin) > l/Vm by causality. We define the time rj V i SC to 
be the moment when the damping scale becomes smaller 
than r]i n , kd(r) V i sc ) = \jr\i n . From that time on, the 
function / decays like a power law, 

f 2 (k, V ) cx k 2 d n + 3 cx / 2 (fc,r7 m )(Wr?r (2 " +3) , 

where a is a positive power describing the growth of 
the viscosity damping scale. Hence, the source term of 
Eq. ( |l6| ) starts to decay faster than 1/a 2 , and additional 
gravity wave production after rj V i SC is sub-dominant. We 
neglect it in our attempt to derive an upper limit for 
primordial magnetic fields. For n > —3/2, the grav- 
ity wave solution given in Appendix Eq. flC5| ) is then 
simply modified by ~log(x m ) — > log(x visc /x m ), since 
the integral of the gravity wave source term only has 
to extend from Xi n to x V i SC . Taking also into account 
that up to rjvisc the cutoff scale is k c [rf) = l/r)i n , hence 
/ 2 (fc, rj) oc k 2 c n+3 = l/?7-™ +3 , we obtain 



dfl G (k) 
dlo^fc) 



( dn B (k) \ 2 {] v 

^dlog(fe) J y^'hn) 



-3-2n 



rad 



-24 log (%viscf •Ein) i 



4 



for n > -3/2 

8(n + 3) 2 log 2 (7 ?1 , 
-3/2 . 



(.Win) 
^rad 
for Tl > 



zf Win) 



(25) 
(26) 



In Appendix we estimate for the two examples of in- 
flation, T in ~ 10 15 GeV, r]i n ~ 8 X 10~ 9 sec and the elec- 
troweak phase transition, T ew ~ 200GeV, ?7j n = ^ eu , ~ 
4 x 10 4 sec, 

Vvisc/Vin — 10 9 for inflation r/i„ = 8x 10 _9 sec, 
Vvisc/Vew ~ 3000 for ew. trans. ?7 et0 ~4x 10 4 sec . 

Up to logarithms, the final formula for gravity wave pro- 
duction is nearly the same for all values of the spectral 
index (cf. Eqs. (§|) and @). 

In these formulas back-reaction, namely the decrease 
of magnetic field energy due to the emission of grav- 
ity waves, is not included. Therefore Eqs. (|23| , |2~4| ) and 

( p5| , p6[ ) are reasonable approximations only if Qg ~ Qb- 
In the opposite case, which is realized whenever 



= { TO) fOT « < - 3 / 2 

" I 8(n + 3)Mog"(U e /^) fOT n > - 3 / 2 



3.3xl0" 6 /i; 



for n < -3/2 



(n+3) 
5xW- e h~ 



(27) 



the magnetic field energy is fully converted into grav- 
ity waves. Note, however, that the value Q,BG( n ) is in 
general not very much smaller than f2 ra d, which is an 
intrinsic limit on Qg for our perturbative approach. 

In Fig. [l| the values Qg and ^B(Vin) as functions of the 
spectral index are shown for two different choices of the 
creation time for the primordial magnetic field: the elec- 
troweak transition, r)i n — r\ ew ~4x 10 4 sec and inflation 
with r] in ~ 8 x 10 _9 sec, for a magnetic field amplitude 
B\ = 10~ 20 Gauss. They are compared with the nu- 
cleosynthesis limit, which comes from the fact that an 
additional energy density may not change the expansion 
law during nucleosynthesis in a way which would spoil 
the agreement of the calculated Helium abundance with 
the observed value. The maximum allowed additional 
energy density is given by 



(28) 
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FIG. 1. We show fic^o an d ^s( 7 7in)^o as functions of the 
spectral index n for two different times of primordial magnetic 
field creation: the electroweak transition ( SJc/io dash-dotted, 
blue and f2s(??m)/io short-dashed, red), and inflation ( Og/i 2 
dotted, blue and f2s(»?in)/io long-dashed, red) for a fiducial 
field strength B\ = 10~ 20 Gauss at A = O.lMpc. The nu- 
cleosynthesis limit, S7ii m /io is also indicated. (The log- terms 
have been neglected.) Clearly, the regimes with Qb > 1 or 
Qg > 1 are not physical and are just shown for illustration. 
We have also shown H,B(Vnuc)h 2 , the magnetic field density 
which is simply cut off at the nucleosynthesis damping scale 
(fat solid line). 

From Fig. [j] we see that fie as calculated above domi- 
nates over flBiVin) for all spectral indices n > —2 in the 
inflationary case and n > for electroweak magnetic field 
production, for an amplitude of B\ = 10~ 20 Gauss. This 
is due to the fact that we have neglected back-reaction 
which leads to a loss of magnetic field energy. Clearly, 
the magnetic field cannot convert more than all its en- 
ergy into gravity waves. However, if our formula for Qq 
leads to Qg > ^B(Vm), it does actually convert most of 
its energy into gravity waves, before it is dissipated by 
plasma viscosity, since gravity wave production happens 
before and at horizon crossing, while viscosity damping 
is active only on scales which are well inside the horizon. 
We can take into account back-reaction by simply setting 
Qg ~ ^B{Vin) when our calculation gives Qq > QB(j]in)- 
We shall use this approximation for Qq m what follows. 

Fig. |l| also shows that, since the value of the magnetic 
field density parameter at which conversion into gravity 
waves is quasi complete is so close to the nucleosynthesis 
limit, ilBG( n )ho ~ 1-12 x 10~ 6 = Qn m h>o, the two curves 
flch-Q and QBiVin)^ cross close to Qn m hQ. This means 
that the gravity wave limit for magnetic fields is very 
close to the limit obtained by setting Qq = ^s(^m)- 

Let us discuss the problem of back-reaction in more de- 



tail. Even if Qg < ^s(?7m), as soon as 



dQg(k) dQ B (k) 
dlog(fc) > dlog(fc) 



5 



for a given scale k 1 , we can no longer neglect back- 
reaction for this scale. The spectrum of fig is 



dfl G (k) 
dbgTfc) 



oc 



k 2n+6 
fc 3 , ' 



for n < -3/2 
for n > -3/2, 



Hence for — 3 < n < 0, the gravity 



Wlllle dlog(fc) 00 fc 

wave spectrum is bluer than the magnetic field spectrum. 
Since there is no infrared cutoff, at sufficiently low values 



■? B ffl and back reac- 



of k we will always have dlo g( fc j < 
tion is unimportant at low k. The value fcn m , below which 
this is the case, can be determined from Eqs. ( plpg| ) and 
©. We find 



fci im A(log 2 (fci im 77 m ))"+ 3 ~ (7^^ 

~ [^^(lO-^Gauss/SA) 2 ]^^, 
for - 3 < n < -3/2 



V2 



(29) 



^lim 7 / 




V8Cl T 



ad 



24fi,„ log 2 (?7„ isc /?/;„)_ 
2 x 10 4 (10" 9 Gauss/B m ) 2 " 



log {Vvisc/Vin) 

for - 3/2 < n< 0, 



(30) 



where 



n A = Bj(/(87r/>c) 



dn B (k) 



and 



fc=l/A 



nx^in/A)' 



n+3 



dbi(fe) 



Bi = Bl(A/r ?in )»+ 3 . 

If fcii m > 1/rjin, e.g. if the square bracket in Eq. ( pQ ) 
is larger than unity, back- reaction is never important. 

For n = the magnetic field and gravity wave energy 
densities have the same spectral index and the condi- 
tion that gravity wave back-reaction becomes important 
is scale independent. In this case it simply reads 



> 



^radVS 



24 log (rivUc/Vu 



(31) 



The situation is different for n > 0. Then the gravity 
wave spectrum is less blue than the magnetic field spec- 
trum and back reaction is always important at sufficiently 
low k, large scales. 

When back reaction is important, it leads to damp- 
ing of the primordial magnetic fields on large scales and 
will actually damp the field down to values for which 
back-reaction is unimportant. This can be seen as fol- 
lows: gravity wave production takes place until Iljj(fe), 
the tensor component of the magnetic field stress tensor, 



vanishes. But then f 2 (k) 
Eq. ph 



which implies according to 



B 2 (q)B 2 (\k- q|) = for all < q < k c . 

For n < the quadratic nature of the coupling of B to 
gravity waves actually damps the magnetic field energy 
at least on all wave numbers q > ky lm /2. 

For n > 0, back-reaction reduces Hij(k) oc 
J d 3 qB 2 (q)B 2 (\k — q|) for small enough values of k. In 
the limit k — > 0, this indicates that back-reaction damps 
the magnetic field on all scales until it becomes unim- 
portant. It is difficult to decide without a detailed calcu- 
lation how the magnetic field spectrum will actually be 
affected, but it seems reasonable to assume that back- 
reaction will alter it until n ~ and the amplitude until 



inequality (31) is violated. We can therefore assume that 
in late time magnetic fields inequality (31) is always vio- 
lated if the magnetic field spectral index is n ~ 0. 

We find this a very important result, which can be 
summarized as follows: Magnetic fields on super-horizon 
scale with a density which is sufficiently close to the ra- 
diation density are strongly damped into gravity waves 
when they enter the horizon. Note also that 'sufficiently 
close' can even mean several orders of magnitude smaller 
since log 2 (kr]i n ) can easily become of order 100 or more. 
Furthermore, primordial magnetic fields produced on su- 
per horizon scales have their spectral index changed by 

gravity wave production to n ~ once they enter the 
horizon. 

During the matter dominated era gravity wave produc- 
tion is somewhat less efficient |Q; and since the scales 
of interest for us are sub-horizon in the matter era we do 
not discuss it here. 



IV. LIMITS AND CONCLUSIONS 

The first limit for primordial magnetic fields produced 
before nucleosynthesis is simply that the energy density 
which they contribute may not change the expansion law 
during nucleosynthesis. As already mentioned, this con- 
dition implies |22] 



£l B (r]nuc)h 2 < 1.12 x 10" 



h 2 

lim^o 



Here we have disregarded the loss of magnetic field energy 
into gravity waves which will, as we shall see, strengthen 
the limit considerably. From Eq. (E2I) we have 



n+3 



4.5/i o z xl0 



B\ (k c (rj nuc )\) 
' 87rp c 2 ^r(^) 

2„m-13( 5 .9xl0 6 )™ 



where we have inserted 

kdijlnuc) 



10- 20 G 



A 



10 13 sec 



n+3 



y/2<7 T tt b p c /(r]3 uc m p ri lad H§) ~ !0 5 /r] n 



G 



~ 6 x 10~ 7 sec^ 1 (for details see Appendix |X| and 
Refs. p9| , p^ , p0[ ). The density parameter Os(ry nuc ) as a 
function of the spectral index n is shown in Fig. [j]. 

Together with the above constraint, this gives already 
an interesting limit on primordial magnetic fields with 
spectral indices n > —2, as shown in Fig. || (solid line). 
For causal mechanisms of seed field production, n > 2, it 
even implies B\ < 10~ 22 Gauss. 



3 

CO 

a 




FIG. 2. We show the nucleosynthesis limit on B\ (solid 
line) as function of the spectral index, n together with the 
limit from gravity waves if the primordial field is produced at 
the electroweak transition (short-dashed) or during inflation 
(long-dashed) for A = 0.1h _1 Mpc~ I0 13 sec. 

Nevertheless, the limit implied from the production 
of gravity waves is more stringent, since the gravity 
waves have been produced at very early times, when 
the magnetic field damping scale was much smaller than 
l/kd(r)nuc) ~ 1.7 x 10 6 sec. The production of gravity 
waves has prevented the magnetic field energy from be- 
ing lost by viscosity damping, since gravity waves do not 
interact with matter in any substantial way. 

Setting Q,q = ^s(?7m) whenever the result of 
Eqs. ( p3j , p6| ) is larger than this limit, which is the simplest 
way to account for back-reaction, the condition 



QgK < 1-12 x l(T fa = Clumh 



lim'tQ 



(32) 

yields the constraint for primordial magnetic fields cre- 
ated at rji n . For spectral indices 



n > -3 



..I 



8f2iii 



-1 



the value for Qq inferred from Eq. ( |26| ) becomes larger 
than Q,B(i]in) at the limiting value On m imposed from nu- 
cleosynthesis (in this approximation we have neglected 
the factor log 2 (j] V i SC / rjin), which can be considerable!). 
Then the magnetic field damping due to gravity wave 



productions is very important. But also for smaller val- 
ues of the spectral index, n > —3, we have flc ~ ^s(?7m) 
for fie ~ fhi m and there is still a considerable amount of 
magnetic field damping due to gravity wave production. 

The results for primordial magnetic fields produced at 
inflation and at the electroweak scale are shown in Fig. ^ 
(dashed lines). As can be seen for the two examples, pri- 
mordial magnetic fields produced before nucleosynthesis 
are very strongly constrained. For all values of the spec- 
tral index, the following expression is a good approxima- 
tion for the limit obtained: 



5 A /l(T 9 Gauss < 700/i x (?7 m /A) (ll+3)/2 AA(n) (33) 

1 . 



where M(n) = J 2^T ' " 



This nucleosynthesis bound becomes stronger for smaller 
cutoff scales, larger k c , according to Eq. (|33| ) it scales like 
(/c c A) - (™ +3 '/ 2 . (Remember that we have set k c = l/?/m.) 

If the seed field is produced during an inflationary 
phase at GUT scale temperatures, where conformal in- 
variance can be broken e.g. by the presence of a dilaton, 
the induced fields must be smaller than B\ ~ 10 _20 Gauss 
for n > —2. If seed fields are produced after inflation, 
their spectrum is constrained by causality. Deviation 
from a power law with n > 2 can only be produced on 
sub- horizon scales, k > l/rji n . Therefore our limit de- 
rived by setting B(k) = on sub-horizon scales, krj in > 1, 
is the most conservative choice consistent with causality. 

Mechanisms which still can produce significant seed 
fields are either 'ordinary' inflation, if the spectral index 
ri ~ —2 or a late inflationary phase at the electroweak 
scale (or even later) where a seed field with n *~ can 
have amplitudes of B\ ~ 10~ 20 Gauss. 

We also have found that magnetic fields which con- 
tribute an energy density close to the nucleosynthesis 
bound, loose a considerable amount (if not all) of their 
energy into gravity waves, which might be detectable. 
In fact, the space born interferometer approved by the 
European Space Agency and NASA, the Large Interfer- 
ometer Space Antenna (LISA) which has its most sen- 
sitive regime where it can detect JIg^o ~ 10 -11 around 
10 _3 Hz~ l/rj wea k 1 22 1 will either detect or rule out all 

magnetic seed fields with spectral index n ~ —0.5 pro- 
duced around or before the electroweak phase transition. 
If LISA does not detect a gravity wave background, the 
constraint analogous to Eq. ([53]) for rj in < 4 x 10 4 sec 
yields 



B\ < 10 Gauss for all indices n > -0.5 

for all mechanisms producing seed fields before or at the 
electroweak phase transition. 

We conclude that, most probably, magnetic seed fields 
have to be produced relatively late, or after nucleosyn- 
thesis to evade the discussed bounds. Our gravity wave 
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bound is not relevant for magnetic fields which are pro- 
duced on sub-horizon scales. But for A ~ O.lMpc to 
enter the horizon, this requires a temperature of creation 
T < IkeV. The only late time mechanism found so far 
which could lead to seed fields is recombination, where 
large scale fields of the order of B ~ 10~ 20 Gauss can 
be induced by magneto-hydrodynamic effects, and the 
difference in the viscosity of electrons and ions p3| , a 
charge separation mechanism. Our work strongly con- 
strains processes of quantum particle production (during 
e.g. an inflationary phase) as origin for the observed 
magnetic fields and favors more conventional processes 
like charge separation in the late universe. 
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APPENDIX A: DAMPING OF MAGNETIC 
FIELDS BY VISCOSITY 

In this appendix we determine the cutoff function 
kd(j]). We use the results found in [p|Jl8| and {jo|. 

We split the magnetic field into a high frequency and a 
low frequency component, separated by the Alfven scale, 
Xa = vat], where the Alfven velocity 



(B 2 



47r(p r +p r ) 



(B\) 



depends on the low frequency component: 
(B w (x)^(x))| Aa , v a ~ 4x I0- 4 x (S A /10- 9 Gauss) 
|L4j . The amplitude of the high frequency component 
then obeys a damped harmonic oscillator equation, with 
damping coefficient, D(r]), depending on time and on the 
mean free path of the diffusing particles giving rise to 
viscosity ]19|| . In the oscillatory regime, we define the 
damping scale at each time rj to be the scale at which 
one e-fold of damping has occurred: JJ' -ydr? = 1. The 
damping term D is given by D — k 2 A co ; / 0,(77) , where 
X co i is the mean free path of the particle species with 
the highest viscosity which is still sufficiently strongly 
coupled to the magnetic field. Long wave modes with 
1/k > va?] are not significantly damped. We now deter- 
mine the damping scale as a function of time. To deter- 
mine whether a given mode with k > kd(rj) is effectively 
damped one has to decide whether it is in the oscillatory 
regime, ujq — kvA > D = k 2 X co i/a(rj) where damping 
really has time to occur or in the 'over-damped' regime 
k < VAa(rj)/2X co i where amplitudes remain approxima- 
tively constant. With va this depends on the magnetic 
field under consideration. 



Let us now determine the damping scale. Before neu- 
trino decoupling at T ~ IMeV corresponding to 7/ ~ 
I0 10 sec, damping is due to both photon and neutrino 
viscosity. The mean free path of photons is 



A 



f 



coZ,7 



a 3 (1.5 x I0 20 sec) , 



where <7t = 6.65 x 10 _25 cm 2 is the cross section of Thom- 
son scattering. For neutrinos, we take into account scat- 
tering with leptons as the principle scattering process 
giving rise to viscosity: 



A 



1 



col.v 



a 5 (7 x 10 48 sec) , 



where a w = G 2 F T 2 is the weak cross section and Gp — 
(293GeV)~ 2 is Fermi's constant. Note that we set Ti = 
c = 1 so that a cross section also can have the units 
GeV- 2 . 

Using the expression for the scale factor given in 
Eq. (^), one finds that photon viscosity dominates un- 
til j] ~ 10 5 sec, leading to 



k d (jj) ~ (2 x I0 10 sec 1/2 )rT 3 / 2 



(Al) 



For j] > I0 5 sec neutrinos viscosity takes over, with cutoff 
function 



k d (v) = (4 x I0 15 sec 3/2 )rr 5/2 



(A2) 



during the oscillatory regime. The comoving wavenum- 
ber k is given here in units of sec -1 . 

After r\ ~ 1 10 sec neutrinos decouple and the dominant 
viscosity is again photon viscosity leading to the cutoff 
function (Al). 



Estimating the viscosity time, namely k ( i(rj V i SC 



f° r inflation, r\i n ~ 10 8 sec and the electroweak 
phase transition, r\i n = rj ew ~ 4 x 10 4 sec, we find from 

the expressions above 1]visc / Tjinl inflation rv ' 

3 x 10 9 and 

Vvisc/Vew ~ 3000. The first result is calculated using pho- 
ton viscosity is just approximative, since we do not know 
the relevant cross sections up to the scale of inflation, 
10 15 GeV, but we certainly expect the value to be very 
large, since interactions are strong and thus viscosity is 
weak. The electroweak result, calculated using the neu- 
trino viscosity, would be quite reliable in the oscillatory 
regime. However, for magnetic fields B < 10 _9 Gauss, 
for which the Alfven velocity is smaller than 10~ 4 , the 
scale r] v isc is still in the over-damped regime. The time at 
which the scale can then effectively be damped depends 
on the value of the magnetic field. In this sense our re- 
sult is only a lower limit, i] V i SC /r) ew ~ 3000. This is not 
very important for our final bounds, where we will even 
set \ogr} V isc/f]in ~ 1, in order to obtain results which are 
independent of the time of magnetic field creation. 

As an example we also determine the damping scale at 
nucleosynthesis, T ~ O.lMeV, z nuc ~ 4 x 10 s which we 
need in Section 4. Setting D-q/2 — 1, we obtain 
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kd(Vnuc) = [2a(ri nuc )a T n e {r] nuc )/r] nuc \ 1/2 



(A3) (ir^k^IT^k'^)) - V l \ b V lm cd {T a \\i,ri)T^ ^ 



Using n e = p c flb / (m p a 3 ) , where m p is the proton mass, 
as well as our expression for the scale factor one obtains 

kdiVnuc) 10" 7 sec _1 ~ 10 5 /?7 imc 

This can of course also be obtained by simply using 
r\ nuc ~ 10 11 se c in the above function for photon viscosity 



given in Eq. ( |Al|) . Again, whether or not this scale is 
in the oscillatory regime and can be effectively damped, 
depends on the value of B(kd)- For B(kd) ~ 10 _6 Gauss, 
which satisfies the nucleosynthesis bound, this is largely 
the case, and for magnetic fields of interest to us kd(j] nuc ) 
is the correct damping scale. 

At the end of the radiation dominated era, photons 
decouple and viscosity acts no more. Since gravity wave 
production in the matter dominated regime is not im- 
portant, we do not calculate the cutoff function in this 
regime. 



APPENDIX B: THE GRAVITY WAVE SOURCE 
OF STOCHASTIC MAGNETIC FIELDS 

The Maxwell stress tensor of a magnetic field in real 
space is given by 

T^(x,^)-i-[ J B i (x,r,)i?^(x,, 7 ) 

-is«(x,7 ? )B B (x,»7)B»(x,» 7 )] . 

In Fourier space, using the Fourier transform convention 
adopted in this paper and the scaling of the magnetic 
field with time, we have 

^(k ) ,) = i -i_/A^( q )^(k- q) 



(Bl) 



where we have introduced the factor 1/a 6 to transform 
the present field B l (k) = B l (k, r]o) back to the physi- 
cal field B l (k, rj) = B l (k)/a 3 . LT^ (k, 77) is the transverse 
traceless component of T lJ (k, 77) , which sources gravity 
waves. Here we give the details of the calculation of 
its correlation function, (IP- 7 (k, 77)IT*' m (k'?7)) which we 
use to compute the induced gravity waves. The projec- 
tor onto the component of a vector transverse to k is 
Pij — Sij — kjtj. Consequently P % a Pl projects onto the 
transverse component of a tensor. To obtain the trans- 
verse traceless component we still have to subtract the 
trace. Hence defining the projector 



' ab 



P* pi pij p , 



(B2) 

To simplify the calculation, we note that up to a trace, 
which anyway vanishes in the projection (|B2|) , T ab (k,r/) 
is just given by 

A " 6(M) S to{l)*a* / AS a (q)B b (k~ q) . (B3) 
We therefore can write 

<ir^,77)ir' m (k',77)) =V l3 ab V lm cd (A ab (k, v )A* cd (k', v )) . 

(B4) 

To co mp ute the two point correlator of A, we use expres- 
sion (B3) and the assumption that the random magnetic 
field be Gaussian, so that we can apply Wick's theorem. 
In other words, products of four magnetic fields can be 
reduced by 

(B*(k)B*\q)B n (s)B* m (p)) = 

(B i (k)B^(q)){B n ( S )B* m (p)) + 
(B*(k)B n (s))(B*i(c l )B* m (p)) + 
(B i (k)B* m (p))(B n (s)B^(c i )) . 



(B5) 



Using also the reality condition, B* a (k) = B a (— k), 
and the two point correlator (||), we obtain 



(A Qb (k, ?7)A* cd (k',ry)) 



d 3 qd 3 p[5(k)S(k') x 



4(2tt) 

B 2 {q)B 2 (-p){S ab - q a g b )(S cd - p c p d ) + 

- p)S(k - q k' + p)B 2 (q)B 2 (\k - q|) x 

(,5 ac - q a q C )(S M ~ (k^q) b (k^q) d ) + 
+<5(q - k' + p)5(k - q p)B 2 {q)B 2 {\k - q|) x 
(S ad - q a q d )(S bc - (k - q) 6 (k - q) c )]. (B6) 

The first term only contributes an uninteresting constant 
and can be disregarded. For the remaining two terms 
integration over d 3 p eliminates one of the two (5-functions 
and leads to 

(A^k^A^k',^)) = 

a- 12 r 

6(k-k')^-^ j d 3 qB 2 (q)B 2 (\k- q \)x 



(^-gTX^-tk-q) ( k -q) ) + 



(6 ad - q a q d )(S bc - (k-q) (k - q) C ) 



(B7) 



we have 



Clearly, the correlator of A and thus also the one of n is 
symmetric in k and k' and hence also under the exchange 
of the first and the second pair of indices. In addition it 
is symmetric in the first and the second as well as in the 
third and the fourth index. The most general isotropic 
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transverse traceless fourth rank tensor which obeys these 
symmetries has the tensorial structure 

M ijlm {\t) = 6 a 8 jm + S" n S jl - S ij S lm + k- 2 (5 ij k l k m + 
S lm k l k j ~ 5 im k j k l ~ S'ttfk" 1 - S jl k l k m 
-5 3m k l k l ) + k- A k l Vk l k m . (B8) 

We could not find a straight forward derivation of this 
result in a textbook on multi-linear algebra where it ac- 
tually belongs, but it can be found, e.g. in [ p4[ . 
We can hence set 

(IF (k, V )Tl* lm (k', rj)) = f(k, r;) 2 /a 12 A^ y7m (5(k - k') 
with 

(n tf (k,T 7 )ir«(k / ,rO) = ^/(fc,ry) 2 <5(k-k'), (B9) 

To determine the correlator of II it is therefore sufficient 
to calculate its trace. With V id abV^ d = V a ujV^ d = Vabcd, 
(for the last identity we simply use that projectors are 
idem-potent), we have 

(^■(k, 77)11^ (k', 77)) =P abcd (A Qb (k,77)A, crf (k',77)) . 

(BIO) 

A somewhat tedious but straight forward computation 
gives 

p abcd^ § ac _ qa^bd _ (k^q) 6 (k^q) d ) 

+ (S ad - q a q d )(S bc - (k^q) 6 (k^q) c )] = 
1 + (k-(k^q)) 2 + (k-q) 2 + (k.q) 2 (k-(k^q)) 2 . (Bll) 

Setting 7 = k • q and /3 = k • (k — q) , and using the fact 
that the second term transforms into the third one under 
the transformation q — > k — q, we finally obtain 



(n ij (k,r 7 )n*«(k',j?)) = 



4(2tt) ; 



:<5(k-k) x 



J d 3 qB 2 (q)B 2 (\k - q|)(l + 2 7 2 + 7 2 /3 2 ) , (B12) 
which leads to the result for f(k) given in Eq. ([10). 



APPENDIX C: GRAVITATIONAL WAVE 
PRODUCTION 

The equation for gravity wave production due to tensor 
type anisotropic stresses is 



hij + 2-htj + k A h l3 = 8nGUi 



(CI) 



For each mode we therefore have an equation of the form 



h + 2-h + k 2 h = s(k,r]) , 
a 



(C2) 



^r-f(k,T)). The function / only depends 
-3/2 via the damping cutoff kd{rj). In 



where s(k, rj) = 
on 77 for n > 

terms of the dimensionless variable x — kr\ equation (C2) 
reduces to 



h" + 2-ti + h = s(k,r})k- 
x 



(C3) 



where a = 1 in the radiation dominated era, and a = 2 
in the matter dominated era. The homogeneous solu- 



tions of Eq. (C3) are the spherical Bessel functions jo , yo 
in the radiation dominated era, and ji/x, yi/x in the 
matter dominated era respectively. We assume that the 
magnetic fields were created in the radiation dominated 
epoch, at redshift Zi n . Using the Wronskian method, the 
general solution of Eq. (C3) which vanishes at Zi n is given 
by 



h(x) = a(x)gi(x) +c 2 (x)g 2 (x) 



(C4) 



where g\ , gi are the above mentioned homogeneous solu- 
tions and 

d(x) = -k- 2 f s(x')g 2 (x')/W(x')dx' 
c 2 (x) = k- 2 f s{x') gi {x')/W(x')dx' , 



W = gig' 2 — g[g 2 is the Wronskian determinant of the ho- 
mogeneous solution. Inside the horizon the homogeneous 
solutions (?i and 172 begin to oscillate. The contribution 
to the integral from times where the scale under consid- 
eration is sub-horizon is hence negligible. Furthermore, 
since the gravity wave energy is growing with wave num- 
ber (it is proportional to k 3 f 2 ), our limit will come from 
large wave numbers, small scales, which ent er t he horizon 
before decoupling. Let us thus solve Eq. (C3) explicitly 
in the radiation dominated regime, rj < rj eq , for a wave 
number which enters the horizon in the radiation era, 
krj eq > 1, and in the case where / is not time depen- 
dent (n < —3/2). We first notice that the Wronskian 
W{jo,yo) — 1/x 2 . Using the radiation approximation 
of Eq. (g) for the scale factor, a = iJo^V^rad we have 

k- 2 s(x) _ 8irGf(k) 



W{x) 



Since y$ diverges at small x the term C\ clearly domi- 
nates. After horizon crossing we have 



sin x 



h(x) ~ cx{l)j (x) = ci(l) 

X 

Performing the integral Ci(l), we find 
, . . 87rG/(fc) sin a; 

h \ X > ~ TT2C, k>g(£m) 



(C5) 



for x > 1 and 77 < rj eq ~ %/^rad/^o- We have com- 
pared this formula with the numerical solution and, as 
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expected, found that it is a very reasonable approxima- 
tion (within less than 10% of the numerical result). 

After horizon crossing, the gravity waves thus prop- 
agate freely, and their energy just scales like radiation 
energy, so that for kr] eq > 1, using Eq. ( |l9|) 



dtl G (k) d PG {k) 



k 3 h 2 



d\og(k) p rad dlog(fc)"" iad a 2 p rad (27r) 6 G 
During the radiation era, on sub-horizon scales 

; 87rG/(fc) „ 3 



9. 



rati 



(C6) 



so that 

dn G {k) 



log(x in ) cos(x) and a p Tad 



BttG 



4fc 3 /(fc) 2 (8ttG) 2 log 2 (x m ) cos 2 (x) 

H 4 r! rad 3(27r) 5 
12fc 3 /(fc) 2 log 2 (x m ) 
/ 9 2 fi rad (27r) 5 



(C7) 



Since the ratio between the gravity wave energy density 
and the radiation energy density is time independent, this 
formula is valid also in the matter era. p c = 3Hq / (8ttG) 
denotes the critical density today. 
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